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r*j ! Abstract 
O-i 

We compute the scale dependence of /nl for models of multi-field inflation, allowing 
for an arbitrary field space metric. We show that, in addition to multi-field effects and self 
interactions, the curved field space metric provides another source of scale dependence, 
which arises from the field-space Riemann curvature tensor and its derivatives. The scale 
dependence may be detectable within the near future if the amplitude of /nl is not too 
far from the current observational bounds. 
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1 Introduction 



Inflation [lj, the dominant paradigm of the early universe, successfully predicts the statistics 
of the primordial perturbations, within the observational limits seen today. However many 
fundamental questions remain, about the form of the inflation Lagrangian, how many fields 
were involved, etc. The upcoming large gains in sensitivity mean that there is a realistic 
chance to provide definite answers to some of them. In particular, non-Gaussianity contains far 
more information than the power spectrum. Non-Gaussianity is most commonly parametrised 
through the bispectrum and for a particular shape (e.g. local or equilateral) all information 
is then contained in a single amplitude, /nl [2]- However, just as is the case for the power 
spectrum, it is rather natural for /nl to be scale dependent. This scale dependence is not 
strongly constrained and therefore could even be much stronger. Provided that the fiducial 
value of /nl is large enough, there is a realistic chance for the Planck satellite to simultaneously 
provide the first measurement of both /nl and its scale dependence parametrised by rif NL in 
the near future. 

One of the aims of early universe physics is to test fundamental theories at far higher 
energy scales than can ever be reached by terrestrial experiment. High energy theories such as 
supersymmetry or string theory typically predict that there are multiple scalar fields, and that 
the kinetic term involves a non-trivial field metric determined by e.g. Kahler potential [3]. This 
field metric affects observable parameters but this topic has not been very extensively studied. 
Most previous work was restricted to the power spectrum, i.e. to linear perturbations. In this 
case it was shown that the field metric gives rise to a new term in the spectral index which 
may easily be as large as the usual slow-roll terms [H El E] ■ Only very recently the first detailed 
study of the bispectrum has been made, following the direction of [7J, in the general case by 
Elliston et al. [8]. They provided covariant formulae for the bispectrum of the field fluctuations 
at horizon crossing and extended the SN formalism [H El El E] to provide a generalisation of 
the well-known SN result for the trivial field space metric [10J to a curved one. There have 
been previous studies of non-Gaussianity from a curved field space metric, but none so general 
or providing explicitly covariant results, see for example [IT] . 

Over the last few years, it has become clear that the scale dependence of local /nl is also 
a sensitive probe of early universe physics, and that observations are sensitive to it. Often the 
scale dependence, n/ NL , is of the same magnitude as the spectral index nn — 1, and in some 
models it may be much larger, for example the self- interacting curvaton scenario [12J. Non- 
Gaussianity may also change rapidly over a short range of scales from being zero on large scales 
to large on smaller scales [13]. Scale dependence of equilateral /nl was first considered in [TJ] 
and for the local model in [T3] (for the specific case of two-field hybrid inflation), with a more 
general formalism developed in [THJ [T7|. Numerous studies of the scale dependence have been 
made for a large variety of models, for an incomplete selection see [18] . A general lesson is: any 
detection of primordial local non-Gaussianity is very valuable since it can rule out all single 
field models which predict a value of /nl in the squeezed limit proportional to the spectral 
index [Tj5] (see however [20]), but for most models in which / NL can be large its amplitude can 
be tuned. Therefore a further observable, either the trispectrum or the scale dependence of 
/nl, will still be required to help discriminate between the large number of models. In this 
paper we will always focus on the local model of non-Gaussianity. 

Sefussati et al. have forecasted that for (local) non-Gaussianity, Planck could reach a 
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sensitivity of a nf ~ 0.1 for a fiducial value of /nl = 50 |2lj. This is not much larger than the 
currently measured central value for the deviation of the spectral index from scale invariance, 
n-ji — 1 ~ —0.04 [22], and shows that there is a real possibility that Planck is able to measure 
both /nl and rif NL together, provided that the fiducial values are large enough. Numerous other 
forecasts have been made especially for surveys including large scale structure (LSS) data [23] 
and of course the tightest constraints may be expected to come from combining the cosmic 
microwave background and LSS data [24]. This is a topical field and very recently the first real 
constraints on n^ NL have been made in [25] (see also [26]). 

This article is structured as follows. In Section [2] we study the evolution of the perturbations 
and rederive the spectral index of the power spectrum as a warm up. In Section [3] we derive a 
general formula for the scale-dependence of /nl, the main result of our paper, and we check our 
result by showing that it reduces to the known result in the case of a trivial field space metric. 
Finally we conclude in Section HJ 

2 Evolution of field fluctuations 

We study a period of inflation driven by an arbitrary number of scalar fields, labeled with 
indices a, b, c, ■ ■ • , with Lagrangian density 



where is the usual space-time metric, and 7^ is the field metric, which is a function of the 
field values. In the case of a trivial field space metric, 7 a & = 5 a b. The potential V is also an 
arbitrary function of field values, except that we assume that it is sufficiently flat so that we 
may use the slow-roll approximation around the time of horizon crossing. 

We take ti to be a fixed, pivot time shortly after all the modes of observational interest have 
crossed the horizon [I], while to < U is the horizon crossing time of a certain mode k = (aH)o, 
and hence k dependent. Hence for each mode we find the number of e- folds between t and t{, 



This situation is depicted in Figure [TJ 

An important point to notice, in order to describe the field fluctuation on the flat slices 
S(j) a = <f) a (t, x) — covariantly, is that in general it is coordinate dependent and thus is not 
covariant. To keep covariance, we realise that (fi a and 0q are uniquely connected by a geodesic 
parametrised by A with respect to the field space metric 7^. Then, 8(j> a can be expressed, with 
(t> a \\=o being identified with the background field, in terms of the vector Q a = d(j) a /d\\\ = o which 
lives in the tangent space at </>q. That is, 



Thus any tensor quantity written in terms of Q a is manifestly covariant by construction [TJ. 
Then, on very large scales where the space-time metric is that of unperturbed Friedmann- 
Robertson- Walker one, using the e-folds N as the time variable the equation of motion for Q a 



C = ~^ lab g^d fi rdu<p b -V, 



(1) 




(2) 



^-^^5r = Q a - l -T a bc Q b Q c + --- . 



(3) 
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Figure 1: Schematic display of how we use the 5N formalism. iVo is the horizon crossing time 
of a mode, while iVj is the fixed "initial" time, which we choose to be shortly after all modes 
of observational interest have crossed the horizon and hence is independent of k. ANk is the 
number of e-foldings between these two times, which obviously depends on k. The final time, 
which we assume to be after all isocurvature modes have decayed and the curvature perturbation 
is conserved, is denoted by Nf. The two early times, N and iV are defined on spatially flat 
hypersurfaces, while the final surface at Nf is comoving [H [27] . 



is® 



D N Q a = w\Q b + \w\ c Q h Q c + 



Wab = U(a;b) + 



Rc(ab)d 0o K 



H H 



W a (bc) = U(a;bc) + ~ 



H H 



D r r^O A p ^0 

ft(a\de\b;c) ~f7~F7 ~ ^^a{bc)d,— 



H 



3H 2 



(4) 
(5) 
(6) 
(7) 



where is a covariant derivative with respect to N, R a b c d,V b = V a - C d — V a -d c and a semicolon 
denotes a covariant derivative with respect to 7^. The use of the parentheses around field 
indices implies symmetrisation, and any terms between vertical bars are excluded from the 
symmetrisation. This implies that 



Q a (N { = N + AN) = Q a (N ) + AN k ( w\Q b + -w\ c Q b Q c + ■■•) + 



Here as always in this text, we will work to leading order in AiV, as this is all we need in 
order to calculate the scale dependence of /nl- Note that AiV itself is order of a few (i.e. of 
order the range of scales which we can observe, which corresponds to 5 - 10 e-foldings), but all 
(AiV) 2 terms will be multiplied by a second order in slow-roll term, unless n/ NL is not small, 
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in which case our formalism breaks down. This is analogous to the situation with the power 
spectrum [5j |6] , one may only treat the spectral index as being small and slowly varying if its 
running is further suppressed. 

From the covariantised version of the 5N formalism |8J, the comoving curvature perturbation 
71 is given by 

TZ k (t f ) = 5N = N a {ti,t f )Q a k {U) + -Nabfatf) [Q a (ti)*Q b (t t )] k + • • • , (9) 

where star denotes a convolution, and the final time tf is assumed to be after all the isocurvature 
perturbations have decayed and 1Z is conserved, so we will drop it from future expressions. In 
the case of a trivial field space metric, the 5N coefficients follow from a direct Taylor series 
expansion, N a = dN/d<f) a and so on [TU]. We will not consider the case of isocurvature modes 
which persist until today in this paper. Note that based on our definition of t iy the derivatives 
of ./V are scale independent. See Appendix B of |T7] for more details, and alternative ways to 
use the SN formalism. 

Now with the fc-dependent factor manifest in ([8]) that comes from the evolution from the 
horizon crossing, which is different for each mode k, to the initial time for the SN formalism, 
we first recapitulate the calculation of the power spectrum Vr and the spectral index n-ji- For 
the power spectrum it is sufficient to work to linear order in Q a (in which case it reduces to 
S<p a ). It is given by 

(K k (t f )K q (t f )) =(27r) 3 <^(fc+ q) 2 -^V n {k) = N a (U)N h (U) (QJfe)Qgfe)) - (10) 
Here, using (jEJ) we can find 

(Ql(U)Q b g (U)) = (Q a k (to)Q b q (to)) + 2AN k w a c (Q b k (t )Q c q (to)) , (11) 
where close to horizon crossing 

(Q a k Q b q ) = ^ {3) (k+q){l ab + e ab ) , (12) 

with e ab being first order in slow-roll and slowly varying jU [5] , and can be found to arbitrary 
higher order in slow- roll using the Green's function approach [6, 28J. Hence its derivative with 
respect to log k is second order in slow-roll. When taking the derivative of (j!2p we need to take 
into account the scale dependence of H 2 (t ) which gives rise to — 2e in the spectral index, but 
the covariant derivative of the metric is zero by definition, so this does not give rise to any scale 
dependence. Similarly for calculating the scale dependence, we may raise and lower indices 
using either 7 ab (tj) or 7 afe (£o) when working to the same level of precision. Finally notice that 
working to the same order, it is only necessary to specify the time dependence of background 
terms, and not those which are multiplied by AiV. 

Hence we may straightforwardly calculate the spectral index from the above expressions, 
and the result is 

^- 1 = ^F = - 2e - 2 nvy^-' (13) 

which agrees with the known result jl]. 
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Notice that this result, together with the observational constraint on the spectral index 
does imply that at least the particular combinations of R a bcd which appear through w ab in the 
expression for the spectral index must be small, barring a chance cancellation between this 
and another term. This provides some justification for our assumption that terms involving 
derivatives of the field metric should be slow-roll suppressed. 



3 Bispectrum and the running of /nl 

The bispectrum of 1Z is defined by 

{n kl {t f )TZ k2 {t f )TZ k ,{t f )) = (27r) 3 5^(k 1 + k 2 + h)B n (k 1 ,k 2 ,k 3 ). (14) 

We will make the usual assumption that the bispectrum due to the non-Gaussianity of the fields 
at horizon crossing [S] does not give rise to observable values of /nl, which may be the case 
for DBI fields and non-Bunch-Davies vacuum. This is generally the case although there 
can be exceptions if the third derivative of the potential is large [30]. There is no analogous 
proof in the case of a non-trivial field metric, but since we are assuming the Riemann tensor is 
slow-roll suppressed this assumption is likely to remain valid, at least in the vast majority of 
cases. 

From ([H]), the three point function consists of two terms, 

(1l kl (t f )n k2 (t f )n k ,(t f )) =N a N b N c (Q a kl Q b k2 Q c ks ) + \ {N ab N c N d ( [Q a * Q b ] fc Q^) + 2 perm} 

(15) 

In the above, all terms on the right hand side (RHS) should be evaluated at £j. We first deal 
with the first term on the RHS of ( fTBT) . Using (jHJ) and ( fP2l) . we can easily find 

N a (U)N b (U)N e (U) {QlMQlMQlM) 

= N a {U)N b {U)N c {U) !^AN kl w a de ( [Q d (t Q )*Q e (t )] ki Qi(t )Q c k3 (t )) + 2 permj 

= (2n) 3 5^(h + k2 + h)N a (t i )N b (t i )N c (t i )^^w abc (k s AN kl + 2 perm) . (16) 



For the second term of the RHS of ( |T5l) . it is useful to use (PHD so it follows that 

^N ab (U)N c (t t )N d (U) ([Q a (U) *Q b (U)] ki Ql 2 {U)Q d ki {t 

= (2n) 3 5^(k 1 + k2 + k s )N ab (t i )N c (t i )N d (t i )?y& ( 7 ac 7 M + 2AN kl w ac 1 bd + 2AN k2l ac w bd ) . 

(17) 

Putting both terms together, and specialising to an equilateral triangle, which is the only 
case we need to consider in order to calculate rif NL (for the justification see [H]), we find that 
/ NL is given by 



6 N ab N a N b 
5 /NL (iV c iV<f 



N d N e N f w def N de N d N f w ef N d N e w de 



AN k - + 4 — - 4 



N gh N9N h N gh N9N h N g N9 
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(18) 



Therefore iif NL is simply the term multiplying AN^ in (|18|) . i.e. 

_ 1 -D/nl _ N a N b N c w abc , ab (N a N b N ac N b N c \ 

This general formula for the observable n/ NL is the main result of this paper. The first term in 
( 1191) is due to the non-linearity of the field fluctuations which are generated between the time 
to and ti, i.e. the correlator (5<f) 3 (ti)). Meanwhile, the remaining term of ( 1T91) arises from the 
effect of multiple fields contributing to 1Z. Note that in the curved field space case, the effects 
of 7° 6 enter into both terms, which makes their separation less clear in this case than in (1201) 
below. 

As can be read from the coefficients w ab and w abc , naively rif NL = 0(e). However since the 
derivative term does not appear in the formula for the spectral index, it could be larger and its 
value is model dependent. It is possible to have a large value of n/ NL [31], which should be a 
sharp prediction of such a model. It is estimated that for a fiducial value of /nl = 50, Planck 
can reach a sensitivity of An/ NL = 0.05 depending on the sky coverage |21j . 

As a check on our result (fT9"|) . and to gain intuition, we demonstrate that our result reduces 
to the known expression in the case of a trivial field space metric. From [T7] , 



n W = _FlN°Njc +l(2 ^_ } (N^ _ W , _ (20) 
/NL N d N e N de V V a 6 lab! \N d N d N d N de N e ' ' 



where we have defined e a = (V^/V) 2 /2, r\ ab = V ab /V and 

F abl = v^2 (^-4:y/e a e b e c + r] aby /T c + r] bc ^T a + r] cay /T b 



1 V. 



abc 



V2SIPJ 

Using w ab = -(V a /V) tb = 2^e a e d - r] ad and w abc = it follows that (EE]) reduces to ([20]). 



(21) 



4 Conclusions 

We have calculated a general formula for the scale dependence of /nl, for the first time allowing 
for a curved field space metric. The non-trivial field space metric, which appears in the kinetic 
term of the scalar field Lagrangian, gives rise to new terms in n/ NL depending on both the 
Riemann curvature tensor and its derivative with respect to the fields. The derivative terms 
do not appear in the analogous formula for the spectral index and hence could be significantly 
larger, since they are not constrained by observations showing that the spectral index is close 
to scale invariant. 

Our work is motivated by the goal of connecting fundamental, high energy theories to the 
rapidly improving cosmological observations. Fundamental theories often predict the existence 
of large numbers of scalar fields, and of a curved field space metric. In 2013 the first relevant data 
from the Planck satellite will be released, and it could lead to a discovery of non-Gaussianity. 
In this case, constraints on not only its amplitude, but also its scale dependence will go a long 
way towards discriminating between the many models of the early universe, especially when 
combined with other measurements such as the spectral index and tensor-to-scalar ratio. 
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